Transmission gap, Bragg-like reflection, and Goos-Hanchen shifts near the Dirac point 
inside a negative-zero-positive index metamaterial slab 
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Motivated by the realization of the Dirac point (DP) with a double-cone structure for optical field 
in the negative-zero-positive index metamaterial (NZPIM), the reflection, transmission, and Goos- 
Hanchen (GH) shifts inside the NZPIM slab are investigated. Due to the linear Dirac dispersion, 
the transmission as the function of the frequency has a gap, thus the correspond reflection has a 
frequency or wavelength window for the perfect reflection, which is similar to the Bragg reflection 
in the one- dimensional photonic crystals. Near the DP, the associated GH shifts in the transmission 
and reflection can be changed from positive to negative with increasing the wavelength. These 
negative and positive shifts can also be enhanced by transmission resonances, when the frequency 
is far from that at the DP. All these phenomena will lead to some potential applications in the 
integrated optics and optical devices. 
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I. INTRODUCTION 

It is well known that a light beam totally reflected 
from an interface between two dielectric media undergoes 
lateral shift from the position predicted by geometrical 
optics [1]. This phenomenon was referred to as the Goos- 
Hanchen (GH) effect 3 

and was theoretically explained 
firstly by Artmann in 1948 [3]. Up till now, the investiga- 
tions of the GH shifts have been extended to frustrated 
total internal reflection (FTIR) 0, H, [1] , attenuated total 
reflection (ATR) 0,11, partial reflection H, El, [11 , 
and other areas of physics [2| , such as quantum mechan- 
ics [l^, acoustics II4I ■ neutron physics Q^, spintronics 
[T^ . atom optics [17|| and graphene [l^ . 

Graphene has become a subject of intense interest 
[T9I [20j since the graphitic sheet of one-atom thickness 
has been experimentally realized by A. K. Geim et al. in 
2004 ^21*1. The valence electron dynamics in such a truly 
two-dimensional (2D) material is governed by a massless 
Dirac equation. So g raphene exhibits many unique elec- 
tronic properties [iSII, including Klein tunneling 22]. On 
the other hand, the Dirac point (DP) in photonic crys- 
tals (PCs) for the Bloch states l25j, l26| is found 
from the similarity of the photonic bands of the 2D PCs 
with the electronic bands of solids. Several novel opti- 
cal transport properties near the DP have been shown 
in [H, m, , such as conical diffraction [23| , a "pseu- 
dodiffusive" scaling [2^, and the photon's Zitterbewe- 
gung [2^ . Very recently, Wang et. al. [13, 1285 realized 
the DP with a double-cone structure for optical field in 
the negative-zero-positive index metamaterial (NZPIM), 



and further the pseudodifFusive property [27| and Zitter- 
bewegung effect [l^ near the DP inside such optically 
homogenous media. 

The main purpose of this paper is to investigate the 
transmission gap, Bragg-like reflection, and GH shifts 
near the DP inside a NZPIM slab. Due to the linear 
Dirac dispersion, the transmission has the frequency or 
wavelength stopping-band, thus the corresponding reflec- 
tion has a frequency or wavelength window for perfect 
reflection, which is analogous to Bragg-like reflection in 
monolayer graphene barrier [29j. This so-called Bragg- 
like reflection in such a simple NZPIM slab is quite dif- 
ferent from that in the ID PCs (for instance, a stack 
of Bragg mirrors), resulting from the destructive and 
constructive interferences. More interestingly, the asso- 
ciated GH shifts in the reflection and transmission can 
be changed from positive to negative with increasing the 
wavelength near the DP. Also these negative and positive 
shifts can be enhanced by transmission resonances, when 
the frequency is far from that at the DP. All these phe- 
nomena will lead to some potential applications in the in- 
tegrated optics and optical devices, such as frequency or 
wavelength filters and frequency-dependent spatial mod- 
ulator. 



II. MODEL 

For simplicity, we consider a TE-polarized light beam 
with angular frequency uj and incidence angle upon 
the NZPIM slab in the vacuum, as shown in Fig. [U the 
dispersion of the NZPIM has a linear dispersion ^27.] 
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fc(w) = (cj - lud)/vd, 



(1) 



with group velocity vjj = {duj/dk)\i^^i^jj, uju is the fre- 
quency of the DP (corresponding wavelength is A/j = 
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FIG. 1: Schematic diagram of the reflection and transmission 
inside the NZPIM slab configuration, where Sr and St denotes 
the GH shifts for the reflected and transmitted light beams, 
respectively. 



^-nclLOD)-, where two bands touch each other forming a 
double-cone structure. Near the DP, the light transport 
obeys the massless Dirac equation as follows: 







A. 

dx 







vd 



(2) 



where ^ — 



are the eigenfunctions of the 



Ez2{x,y,uj) 

electric fields with the same k{uj). It is noted that the 
condition for realization of the DP in the homogenous op- 
tical medium is the index varying from negative to zero 
and then to positive with frequency [27[, which is called 
as NZPIM. For simplicity, we take the Drude model as 
the parameters for both the relative permittivity and per- 
meabihty of the NZPIM ^28,]: 



(3) 



(4) 



where uj^^ and uj^p 



are the electronic and magnetic 
plasma frequencies, and 7e and 7m are the damping rates 
relating to the absorption of the material. Here we can 
assume 7e = 7e = 7 ^ '^epj'^mp- It is important that 
when ujep — ^mp — and 7 = (no loss), then both 
£i(a;£)) and ^i{ljd) may be zero simultaneously. In this 
case, we find k{'WD) ~ and vd — c/2, where c is the 
light speed in vacuum 28]. In what follows we will dis- 
cuss the reflection, transmission, and the associated GH 
shifts near the DP in NZPIM slab. 



III. REFLECTION AND TRANSMISSION 

In this section, we will firstly investigate the properties 
of the reflection and transmission. Assuming the incident 
plane wave, £'!,"(a;,y) = exp [«(fc^2; -I- fc^y)], where = 



ko cos ( 



fco sin 9o, ko = lu/c is the wave vector 



in vacuum, the reflected and transmitted plane waves 
can be expressed by El^^{x,y) — r exp [i(— /c^-a; -I- /c^y)] 
and El^{x,y) = texp {i[—kx{x — a) + kyy]}, where the 
reflection coefficient r is 
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,(5) 



and the transmission coefficient is t = e^^/g with the 
following complex number. 



ge^'^ = cos kixa - 
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with kix — y kf — ky and ki = {lj~ujd)/vd near the DP. 
It is clear that the wave vector kix depends on the differ- 
ent frequencies lu and parallel wave vector ky , which will 
resulting the unique properties of reflection and trans- 
mission in two cases oi to > ud and uj < ujz). 

Case 1: Lo > ujo- The reflection probability R can be 
given byEq. ©, 



R 
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V 



Mo kix f^l kx \ . 2 1 lr:\ 

1 — sm kixa, (6) 

Ml kx Mo kix. 



and the transmission probability T — 1 — i? is also 
given by T = |tp — Under resonance conditions, 

kixa — Ntt, {N — 0,1,...), the reflection probability 
Rmin reaches the zero and the transmission probability 
Tmax is equal to 1. Otherwise, at the anti-resonances, 
kixa ^ {N + l/2)7r, {N = 0, 1, ...) the reflection proba- 
bility R tends to 
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and the corresponding transmission probability T„ 
equal to 



T _ A / f f^a kix , Ml kx 

^ min — ^/ I 

VMi kx Mo kix 



(7) 

n is 



(8) 



However, we emphasize here that the reflection and 
transmission can be divided into evanescent and propa- 
gating modes, taking the influence of the incidence angle 
9a into account. The propagation of the light beam inside 
the NZPIM slab can be evanescent when > 9c, where 
the critical angle for total reflection can be deflned as 



2 1 



LU J 



(9) 



with the necessary condition lud < lu < 2lu£,. In 
this case, the transmission and reflection probabilities 
damped exponentially in the following form: 



(10) 
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FIG. 2; (Color online) The transmission gap as the function 
of the frequency u, where a = 100 mm, and uid = W x 2n 
GHz. Solid and dashed curves correspond to 6o = 10° and 
00 = 30°. 



FIG. 3: (Color online) The reflection probability R as the 
function of the wavelength A, where 9o = 20°, and other pa- 
rameters are the same as in Fig. (2] Solid and dashed curves 
correspond to a = 100 mm and a = 10 mm. 
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(11) 



where n = [ky — ki]^/^ is the decay constant. As a matter 
of fact, the light beam can transmit though the NZPIM 
slab in propagating mode at any incidence angles, when 
the critical angle 9c is no longer valid for lo > 2ujd- 

Case 2: lo < ud- The reflection and transmission 
probability can be also damped exponentially when the 
incidence angle is larger than the critical angle, 



sm 



2(^-1 



(12) 



with the necessary condition |wz) < < cj/j. On 
the contrary, the reflection and transmission probabili- 
ties will oscillate periodically on the thickness a of the 
slab, as mentioned Case I. In this case, the refractive in- 
dex, defined as ni = — ^£i/ii, should be negative, which 
will lead to the negative GH shifts, as discussed later. 



Based on the mentioned-above properties of 
tion and transmission in these two cases, the 
sion as the function of frequency to has a gap, 
in Fig. [21 where a = 100 mm, and u>d = 
GHz. Sohd and dashed curves correspond to 
and 9o = 30°. Similarly, the transmission 
a stopping-band for the wavelength, since A 
Since kf^ = (w — uid)'^ /vj^ — fc^ < 0, the 
region of the transmission gap in Fig. [2] is 
— kyVD < uj < ujd + kyVD, which leads to 
of transmission gap as follows. 

Aw = 2kyVD- 



the reflec- 
transmis- 
as shown 
: 10 X 27r 

00 = 10° 
also have 
~ 2ttuj/c. 
frequency 
given by 
the width 



(13) 



This means Alu/uj = sin^o with the help of ~ c/2. It 
is further shown that the transmission gap with the cen- 
ter uj = ud becomes narrower with the decrease of the 
incidence angle, and even vanishes at normal incidence. 
This transmission gapjwhich is analogous to that in sin- 
gle graphene barrier 29\, is due to the evanescent waves 
in two cases oi uj > lod and lo < lod- Furthermore, the 
tunable transmission gap can be further understood by 
the dependence of the critical angle on the frequency, lo. 



Fig. [3] indicates the dependence of corresponding re- 
flection probability R on the wavelength A = 27rw/c, 
where Oq — 20°, and other parameters are the same as in 
Fig. [21 Solid and dashed curves correspond to a = 100 
mm and a = 10 mm. It is interesting that the light beam 
can be perfectly reflected by such single NZPIM slab at 
some range of the wavelength. As indicated in Fig. [31 
the wavelength window for perfect reflection will become 
narrower with the increase of the width of slab, caused 
by the decay factor exp {—2Ka) in Eqs. (fTUj) and pT|) . It 
is clearly seen from Fig. [21 that the reflection also has a 
similar frequency window for the perfect reflection, since 
R = 1 — r. These frequency or wavelength passing-band 
in reflected discussed here is similar to but different from 
the Bragg reflection in the ID PCs. This so-called Bragg- 
like reflection discussed here is exactly due to the linear 
Dirac dispersion described by Eq. ([1]), which results in 
the evanescent waves in two cases of w > lod and lu < loo, 
corresponding to the two eigenfunctions of electric fields 
with the same k{Lo). In a word, the Bragg-like reflection 
will provide alternative way to realize the frequency or 
wavelength filters with more design flexibility and minia- 
turization. 
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FIG. 4: (Color online) The GH shifts as the function of the 
wavelength, A, where a = 100 mm, and other parameters 
are the same as in Fig. (2] Solid, dashed and dotted curves 
correspond to Oo = 30°, 9q = 20°, and Oo = 10°. 



IV. GOOS-HANCHEN SHIFTS 

Now, we have a look at the GH shifts in the reflection 
and transmission inside the single NZPIM slab. When a 
well-collimated light beam with the central incidence an- 
gle do is considered, the GH shifts in reflection and trans- 
mission, according to Artman's stationary phase method 
Q, can be defined as 



Sr,t 



dky 



(14) 



where ky = fegsin^, 6 represents the incidence angle of 
the plane wave component under consideration, (jj^ = 
(j) + Tr/2 and 4>r — 4> are the phase shifts of the reflected 
and transmitted light beams, respectively. Clearly, the 
GH shift in transmission is the equal to that in reflec- 
tion inside such symmetric slab configuration, because 
the values of the derivation of the phase shifts with re- 
spect to ky are the same. Fig. 2] shows that the GH shifts 
can be positive and negative, where a = 100 mm, and 
other parameters are the same as in Fig. [21 Solid, dashed 
and dotted curves correspond to Oq = 30°, Oq = 20°, and 
9q = 10°. It is shown that the GH shifts can be positive 
for A < Ad , while they can be negative for A > Xd- More 
interestingly, the GH shifts near the DP can change from 
positive to negative with the increase (decrease) of the 
wavelength (frequency) . 



Fig. [3] also shows the dependence of the GH shifts on the 
incident angle where a — 100 mm, and other param- 
eters are the same as in Fig. [2] Solid, dashed, dotted 
and dot-dashed curves correspond to w = 13 x 27r GHz, 
w = 11 X 27r GHz, w = 9 X 27r GHz, and w = 8 x 27r GHz. 
It is reasonable that the GH shifts can be negative in the 



FIG. 5: (Color online) The GH shifts as the function of inci- 
dence angle 6o, where a = 100 mm, and other parameters are 
the same as in Fig. [2] Solid, dashed, dotted and dot-dashed 
curves correspond to = 13 x 27r GHz, lo — 11 x 27r GHz, 
uj = 9x2tt GHz, and = 8 x 27r GHz. 



case of w < uj]j, where the refractive index ni — — ^£i/ii 
is negative, while the GH shifts are positive in the case 
of w > lod, where the refractive index rii = ^/Sifii is 
positive. In addition, it is also shown that the GH shifts 
near the DP have only the order of wavelength due to 
the evanescent waves. The smallness of the GH shifts 
are similar to those in total reflection or FTIR structure. 
However, when the frequency is far from that at the DP, 
the incidence angle will be less than the critical angle, or 
there is no critical angle as discussed above. Thus, the 
negative and positive GH shifts can also be enhanced by 
the transmission resonances, as shown in Figs. |4] and [5l 
In addition, the GH shifts also depends on the width a of 
the slab. It can be predicted from Ref. ^ that the neg- 
ative and positive GH shifts in the evanescent case will 
saturate to a constant with increasing the slab's width, 
when the incidence angle is larger than the critical an- 
gle. In a word, these negative and positive GH shifts 
are applicable to realize the frequency or wavelength fil- 
ters in spatial domain, and frequency-dependent spatial 
modulator. 



V. CONCLUSION 

In conclusion, we have investigated the transmission 
gap, Bragg-like reflection, and the associated GH shifts 
inside the NZPIM slab. It is found that the transmission 
has a frequency stopping-band, thus the corresponding 
reflection has a frequency or wavelength window for the 
perfect reflection. This so-called Bragg-like reflection, 
resulting from the linear Dirac dispersion of NZPIM, is 
similar to but different from the Bragg reflection in the 
ID PCs. In addition, the GH shifts near the DP can be 
changed from positive to negative with increasing wave- 
length, based on the unique properties of the reflection 
and transmissions. These negative and positive shifts can 



5 



also be enhanced by transmission resonances, when the 
frequency is far from that at the DP, c£n. With the ex- 
perimental realization of the NZPIM [30| , we hope these 
phenomena will lead to some applications in the inte- 
grated optics and optical devices. 
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